METHOD AND APPARATUS FOR INTERFERENCE REJECTION 

leSiod^and apparatus for rejection of noise and interference from a received combination signal by 
estimation of the desired signal having in addition to noise and interference been distorted by a 
5 communication channel through which the modulated transmitted desired signal has been passed. 

DESCRIPTION OF BACKGROUND ART 

In digital commxmication systems a transmitted signal is normally encoded, modulated and 
10 transmitted over a commimication channel. The signal is sent in form of pulses that can be detected 
either 0 and 1 and given combinations of these "bits" form symbols v^th different meanings. One or 
more signal sequences sent in a packet data stream is called a "burst". 

^ The information signal is sent by modulation of a carrier wave. Different modulation methods exist, 
V§ of which the most usual are based on amplitude, frequency and/or phase modulation. 
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It is advantageous to describe different modulation forms by a complex signal representation. The 
complex signal represents the real physical bandpass signal in an equivalent lowpass representation. 
Essential features of the different modulation forms can then be represented geometrically. For 
2® example in one-dimensional modulation, the signal constellation points are located on a straight 
rh line. 

Furthermore, there is a difference between linear modulation (different forms of amplitude 
modulation) and constant envelope modulation (angle-, phase and frequency modulation). 
25 However, almost all modulation methods used in mobile telephony can be interpreted as linear 
modulation with some modifications. In linear modulation, the modulator and the demodulator 
carry out some signal handling functions to send information by modulation of a carrier wave. 

Whatever the physical medium used for transmission of the signal is, the essential feature is that the 
30 chaimel distorts the signal in different ways. The distorted signal is further corrupted by noise and 
interference. At the receiver, the signal is reconstructed and demodulated. The reconstructed signal 
is processed by channel and source decoders to obtain a copy of the transmitted source signal. 



1 



• 



In a receiver with only one antenna the demodulator normally tries to combat distortion caused by 
the communication channel and noise but neglects the effects of interference. In order to derive 
efficient algorithms, the noise is normally considered to be white. If the noise is not white, a noise 
whitening filter can be used but this is not normally done since the color of the noise is unknown. 

5 

Filtering structures can be used to perform interference rejection. Adaptive or non-adaptive 
algorithms can be used to adjust the parameters in these filtering structures so that the desired signal 
can be estimated. After the filtering structures, demodulators demodulate the estimated signal of 
interest. 

10 

There are two classes of methods for choosing the parameter values. Conventional methods rely on 
the knowledge of training sequences embedded in the desired signals, this knowledge being used to 

:?>*^ 

choose an initial setting for the coefficient values, typically accordmg to a least-squares error 
\l criterion. In an adaptation by means of a training sequence, the received data are used to find the 
1^ parameters. Blind methods use one or more properties of the interference or of the desired signals 

hi I 

[|1 instead of training sequences. Conventional methods that use a known training signal are well 
known to those skilled in the art of signal processing and are described in the open literature of 
Q adaptive signal processing and adaptive filtering. 

ZQ The desired signal is estimated by a sequence estimator, of which there are several versions known 
r- i in the art. These estimators make use of algorithms to mathematically calculate the desired signal by 
using filter parameters. 

The most common sequence estimator used in receiver structures is the Maximum Likelihood 
25 Sequence Estimator (MLSE). The MLSE is preferably implemented with the so-called Viterbi 
algorithm. A version of this algorithm is the Soft Output Viterbi Algorithm (SOVA), which 
improves the performance of the subsequent channel decoding of the signal. The Maximum A 
Posteriori (MAP) algorithm gives superior performance, but is computationally more difficult. In 
the Decision Feedback Equalizer (DFE), already decided symbols in the process are used to decide 
30 the current symbol. A compromise between MLSE and DFE in performance is the Decision 
Feedback Sequence Estimator (DFSE), where it is possible to make a trade-off between 
computational complexity and performance. In the case with more than one antenna, some receivers 
can also consider the interference. 
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The main drawback with existing one-antenna solutions is that the algorithms consider only noise. 
Meanwhile, many communication systems are interference limited. The problem is that when the 
properties of the interference do not fit the model of the noise for which the receivers were 
designed, the performance of the receivers is greatly degraded. 

5 

The Interference Rejection Combining (IRC) receiver combines the output from antennas while 
trying to minimize the effect of noise and interference. There is an algorithm in the literature that 
considers the interference as a colored process. In deriving this algorithm, it is assumed that the 
noise can be described as a white sequence passed through an Auto Regressive (AR) filter. 

10 

The following patent documents are presented as prior art solutions in interference cancellation by 
means of different filter structures and for estimation of the desired signal by selection of suitable 
filter parameters for the filter structures. 

= as: 

yi 

l| In WO 98/16021 there is presented prior art solution for extracting a signal of interest from plurality 
of spectrally and temporally overlapping signals containing digital data. The apparatus of this 

43 solution comprises different filters producing time shifted output signals, frequency shifted output 
signals and linear combining means for summing said output signals and to produce an estimate of 

r- r said signal of interest. 

20 

f ? WO 97/1 1544 makes use of a noise-predictive maximum-likelihood (NPML) data detection scheme 
operating on signal samples received via an equalizing filter. 

EP 0782260 A2 presents an equalizer configuration for processing real-valued and complex-valued 
25 signal samples. 

The development of the mobile communication today is towards more and more users, which leads 
to increasing interference. This means that it is of great importance to find better^and better methods 
for interference rejection and cancellation all the time and that every advantage within this area is 
30 usefiil. 
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SUMMARY OF THE INVENTION 



The object of the invention is a better method and apparatus of interference rejection that offers 
improvements in performance compared to prior art methods. 

The methods and apparatus of the invention, which choose the values of the parameters in the 
filtering structures, are improved methods compared to prior art methods in that a significantly 
better performance is obtained. The method works with one or more antennas as well as with and 
Mdthout oversampling. 

The method of the invention is mainly characterized in that the signal is received as a combination 
of noise, interference and the distorted desired signal through one or more antennas. The received 
signal is separated into a real and an imaginary part and the noise and interference component is 
modeled as a filtered process. By utilizing the combined structure of the signal thus obtained, an 
equation for the received signal is formed as a function of the desired signal distorted by the 
communication channel and the noise and interference component. Values for the filter parameters 
can then be selected in the equation, where after the desired signal can be estimated by calculation 
from the equation by means of the filter parameters selected. 

The apparatus is mainly characterized by 

a) means for receiving a signal as a combination of noise, interference and the distorted desired 
signal through one or more antennas, 

b) means for separating the received signal into a real and an imaginary part, 

c) means for modeling the noise and interference component as a filtered process, 

d) means for forming an equation for the received signal as a function of the desired signal 
distorted by the communication channel and the noise and interference component by utilizing 
the signal structure obtained in steps b) and c), 

e) means for selecting of values for the filter parameters in said equation, 

f) means for estimating the desired signal by calculation from said equation by means of said 
filter parameters. 



In the case when the symbol to be transmitted is one-dimensionally modulated in the complex 
plane, the present inventix)n solves the interference problem by utilizing this inherent signed 
structure. 

The desired signal has been one-dimensionally modulated by a sender and thereafter transmitted 
through a complex valued communication channel. 

The received signal is known to be a combination of the transmitted desired signal and noise and 
interference. Thus, the distortion can be taken into consideration when estimating the desired signal. 
The distortion caused by the communication channel can be modeled mathematically. This model 
can include the distortion of the desired signal as well as the addition of the interference and noise. 
Thus, the distortion is taken into consideration by modeling the channel as a linear filter in said 
equation so that the output of the filter is formed as a product of the filter vector and the input signal 
vector, like in the following equation 

x(t) = hd(t)^n(t) 

where the diacritical marks denotes that the entity to which it is attached is complex-valued. 

The received signal is separated into a real and an imaginary part, so that the received signal can be 

expressed as a fimction of the time t in the following form 

x(0 = hd(OH-n(0 



where the channel matrix has the following form 




and 



x(0 = 


Re(x(/)) 


, n(0 = 


Re(«(0) 




Im(i(/)). 




Im(«(/)) 




and where x(t) is the received signal, d is the desired signal and n is the combined noise and 
interference term and the matrix h contains the channel coefficients. The channel coefficients model 
the distortion of the signal through the commimication channel. 

5 The noise can be modeled as a filtered process, for example as an AR-process, in which the noise is 
described in the following way 

K 

10 

where A represents the filter parameters, K the number of the filter parameters, / is the time and e is 
'^V: Si white noise vector, 

yi 
?" 

Also other filter processes can be used in the invention, for example the MA-process, the 
]^ mathematical equation of which is known and can be found in the open literature describing filter 
processes. 

7^'' The chaimel coefficients and the filter parameters in the modeled process for the desired signal and 
4^ the noise are calculated by means of a known training sequence from the equation formed by 
2^ knowledge of the data sent by the signal in the training sequence and by knowledge of the statistical 
properties of the noise component e. There exist mathematical algorithms by which h and A can be 
calculated firom the equation above when x is known (being the received signal) and when the 
statistical properties of e are knovm. Examples of such algorithms include least-squares and 
maximum likelihood methods. 

25 

After having selected the filter parameters and the channel coefficients, the desired signal is 
estimated mathematically by means of an equalization method using same type of algorithms 
mentioned above. 

30 The equalization methods are carried out in a known manner for those skilled in the art and can for 
example be carried out by means of a sequence estimator. In mathematical terms, the desired signal 
can for example be estimated by maximum likelihood for the equation giving the received signal, 
whereby the value for the received signal giving the least error can be obtained by means of a 
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mathematical algorithm. When the Viterbi algorithm is used, the desired signal can be obtained 
without being forced to calculate the received signal for all different values of the desired signal. 
The idea of the mathematical approach is that the received signal is recreated by means of the 
channel coefficients and the selected filter parameters and different values for the desired signal, 
5 where after the recreated value of the received signal is compared with the true value for the 
received signal to obtain error values for the received signal. The desired signal is then estimated by 
selecting the sequence giving the least error for the received signal. 

The method of the invention can be repeated for each received symbol sequence. 

10 

As was mentioned before, there are different known filter process models that can be used in the 
invention to model the noise signal. These processes can be selected among filter processes usually 
'''li used to estimate the desired signal itself. In an MA model, the signal is a weighted simi of delayed 
IH input components, while in an AR process the signal is a weighted sirni of foregoing output values. 
IS The noise is preferably modeled as a white noise through an AR filter but also other filter processes 
can be made use of. An example of another alternative is the MA filter process. 

ri The choice of the number of filter coefficients worth to be taken into consideration, i.e. the order of 
the AR process, can be estimated adaptively by means of known algorithms. One simple order 
2|0 estimation approach is to choose the largest order that gives a significant decrease in the prediction 
^ error. Another way to choose the proper model order is to use each estimated model to estimate a 
number of transmitted symbols and the corresponding error. The error for different model orders 
can be calculated recursively, i.e. the parameters calculated in the current step can utilize the 
parameters calculated in the previous step in order to simplify the computational complexity. 

25 

The method of the invention can be made computationally less complex by prefiltering. The basic 
idea with prefiltering is to decrease the time spread of the channel, thus reducing the complexity of 
the sequence estimator. The prefiltering can be made by knovm methods. 

30 It is possible to gain some additional performance if oversampling is used. Oversampling means 
that the sampling is faster than the bit rate or the symbol rate. By introducing oversampling, more 
channels are obtained as we get more samples and mathematically this fact can be put in the 
equation in the term describing the length of the, channel and the amount of samples. 
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Prior knowledge of the filter structures can also be used to reduce the amount of parameters. It is 
then primarily a question about the physical properties of the channel which are utilized in the 
equation. 

Decision direction is also a well known method of improving the performance of the system. The 
concept is based on the idea that once all or some of the transmitted symbols are estimated, the 
values are used in a renewed estimation of the channel and filter parameters. These estimates are 
then used for renewed estimation of the transmitted symbols. 

The received signal is preferably linearly modulated with a one-dimensional constellation in the 
method of the invention, but if the received signal has been non-linearly modulated, the signal can 
in some cases be de-rotated to obtain a linear modulation.. One example of a modulation method 
that can be de-rotated is MSK (Minimum Shift Keying). Another example is GMSK (Gaussian 
MSK), which is the modulation scheme in GSM. Although GMSK is a non-linear modulation 
scheme and cannot be transformed into a linear one, it has been shown that it can be approximated 
to a linear one with good accuracy. 

It has been observed through extensive simulations that in interference environments, the receiver 
described in the present invention offers substantial improvements in performance, primarily in the 
form of lower Bit Error Rates (BER), which in turn allows operation with much lower Carrier to 
Interference Ratios (CIR). 

SIMULATION RESULTS 

The method of the invention is simulated for different values of the order of the AR process, namely 
K=l and K=2. The invention is also simulated for an embodiment of the invention, wherein the 
order is chosen adaptively. The channel is a typical urban channel described in the GSM 
specification, and it is also the GSM system that is used for communication. White Gaussian noise 
is also added with a Eb/NO value of 25 dB, where Eb is the energy per bit and NO is the noise 
spectral power. 

The results of the simulations are shown in Figure 1, wherein the performance of the invention is 
compared to a conventional receiver. The result is shown as the BER as a function of the carrier to 
interference ratio (CIR). The conventional receiver used is a Viterbi algorithm with five taps. This 
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^ conventional receiver is very similar to what is used in most state of the art solutions. It fulfills the 
demands in the GSM specification with a 2 dB margin. At a BER level of 3% the improvement with 
the method of invention is about 7 dB. It should be noted that this is an extremely large 
improvement. Normally an improvement of a fraction of a dB is considered worthwhile. The 
5 computational complexity when K=2 is about four times the conventional receiver. For the case 
K=l the improvement is less than 2 dB. The adaptive method that automatically chooses the order 
of K performs very well also. 

In the following, some preferred embodiments of the invention will be described by means of 
10 figures and particular mathematical calculation algorithms. One skilled in the art understands that 
the invention can also be performed with other equalization methods and filter processes than those 
presented in the figures and that the maximum likelihood criteria equations can vary. In other 
O words, the details of the invention can vary within the scope of the claims. 

i 

BRIEF DESCRIPTION OF THE FIGURES AND PREFERRED EMBODIMENTS 

yl 

%^ FORMULAS 1 - 23 describe an estimation of the desired signal for the preferred embodiment of 

4^ the invention 

h'~ 

20 The performance of the invention compared to a conventional receiver is displayed in Figure 1 . 

^■f The main characteristics of the present invention are schematically described in Figure 2. 

CP 

DETAILED DESCRIPTION OF THE FIGURES AND PREFERRED EMBODIMENTS 

25 

An example of a detailed mathematical model for the estimation of the desired signal will now be 
presented by means of formulas 1 -23. The mathematical calculations themselves are known by one 
skilled in the art and are used in the invention in which the desired signal is estimated by dividing 
the signal in a real and an imaginary part and by modeling the noise as a filtered process. 

30 

A transmitted signal d(t) signal , which is received in the method of the invention has been linearly 
modulated and passed through a complex communication charmel h. The communication channel 
cause a distortion on the desired signal and as, the received signal is a combination of the distorted 
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desired signal, and the noise and interference component, the complex-valued signed x(t) can be 
formed as in equation 1 . 

The diacritical mark denotes that the corresponding entity is complex-valued and n(t) denotes the 
additive noise and interference. In equation 1, h is a vector of the form according to equation 2, 
where the subscript r and / denotes the real and imaginary part of the coefficients, respectively, and 
L is the length of the communication channel in bit periods. 

The vector d(t) in equation 1 contains the delayed samples of the transmission symbols and can be 
expressed as in equation 3 where T denotes the transpose operator. 

By separating equation 1 into its real and imaginary counterparts we get equation 4, where h can be 
expressed as a matrix according to equation 5, andx(t) and n(t) as in equations 6a and 6b, 
respectively. 

It has to be noted that in equation 4, all entities are real-valued. The model in equation 4 has one 
input and two outputs. If h is not degenerate e.g. is identically equal to zero, d(t) can be estimated 
from both the real and imaginary part of the received signal (the first and the second row of 
equation 4). Hence it is clear that if d(t) is real-valued there is an inherent redvuidancy in the 
received signal. This is the redundancy that will be used to design an interference rejecting receiver 

In this embodiment of the invention, we choose to model the noise and interference as an AR 
process according to equation 7, where eft) is a white noise with covariance matrix expressed as in 
equation 8. 

Equation 7 can be rewritten as equation 9, where N(t) can be expressed as in equation 10 and Ma as 
in equation 1 1 . 

An equation for the received signal will now be formed. Equation 4 can be written in an augmented 
form according to equation 12, where X(t) is defined similarly to equation 10 so that D(t) is written 
as in equation 13 and H as in equation 14. 

By multiplying equation 12 with Ma equation 15 is obtained. 
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If Ma and H would be known, the source signal could be estimated by minimizing the maximum 
likelihood criteria according to equation 16. Different equalization methods exist for this 
estimation and usually this is done by means of a so called sequence estimator, known in the art. 
The MLSE (Maximum Likelihood Sequence Estimator) mentioned is implemented by a 
5 mathematical algorithm, preferably the Viterbi (or SOVA) algorithm. 



In practice Ma and H are not known and hence they must be estimated. If we have a training 
sequence, sent in connection with the desired signal and known by the receiver, i.e. d(t) is known 
for certain time-indices, the criterion in equation 17 can be used to find the parameters in Ma and H 
10 as well as in Q. In general this is a nonlinear minimization which can be numerically cumbersome. 
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Embodiments of the invention with other filter structures and other orders of the process will now 
be described. Above the noise and interference has been described as an AR process. If a definition 
according to equation 1 8 is made by not considering the structure of B, the parameters can instead 
be estimated as in equation 19. These estimates and equation 20 can then be used in the estimation 
of the unknown transmitted symbols. The minimization problem in equation 19 is a simple linear 
least square problem, which is well described in the literature. When the simplification in equation 
18 is introduced, we get a so called ARX model for the received signal. By the definition in 
equation 18, we have fixed the order of B to L+K . In one embodiment of the invention B is set to 

2S9 the order, which is not necessarily fixed to Lh-K, Normedly, we would have M> or = L), but the 

St 3 ^ 
invention is not limited to this case. 

There will now be described an embodiment of the invention wherein an adaptive method is used to 
choose the order of the filter process. Assume for instance that the receiver operates in a pure noise 

25 environment, so that A1=A2=0, or equivalently n(t) = eft)- Then it is desirable to set the order of 
-> the AR model to K=0 in order not to describe the process with too many parameters. In other 
words, we want to avoid an over-parameterization of the problem, which might lead to a worse 
situation than the correct model. A version of the invention in which the model order is adaptive. 
There are many methods used to determine the model order available in the literature. Most of them 

30 rely on some criterion that is based on the magniUide^of^e esjimation-eir^ 

reflect how well the model order fits the data used in the estimation. If we increase the model order, 
the magnitude of the residuals will always decrease. There is needed a test telling if this decrease is 
significant or not. In this embodiment of. the invention, the magnitude of the residuals is given by 
the norm of the estimated of the noise covariance matrix Q. Many order estimation procedures have 
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been described in the literature that apply to the present invention. Maximum likelihood methods, 
like the Final Prediction Error (FPE), Akaike's Information criterion (AIC) and its variants (such as 
the so called BIC and CIC) and Minimum Description length (MDL) are all applicable to the 
problem treated in the invention. Also sub optimal methods exist. These methods rely on the fact 
that the covariance matrix of the received signal has a certain rank structure that can be used to infer 
the order of the system to be identified. A major drawback with the order estimation methods is 
their computational complexity, why more practical methods are implemented. One simple 
algorithm that estimates the order of the model is outlined in the following. The algorithm is given 
for the case M=K+L, but can be generalized to general ARX models. 

ri) set K=0, r-i = infinity ^^r— vA^^^yC-^^ 
^ 2) estimate the parameters of the systems v^th the current model order 

3) evaluate rK= IIQII 

4) If FK/nc-i is significantly smaller than 1 and K<Kinax, increase K with one and go to step 2 

5) Choose the K that gave the latest significant reduction of nc 

The term "significantly" in the algorithm above is defined by a certain threshold. This threshold is a 
trade-off between the performance in noise and interference environments. The smaller the 
threshold, the better the performance in noise. 

A conceptually different method to choose the proper model order is to use each estimated model to 
determine a number of transmitted symbols. The resulting metric describing the estimation error 
from the Viterbi algorithm can then be used instead of the residual. The advantage of this method is 
that it provides a model order validation, which is independent of the data used in calculating the 
model. 



Another way of having adaptive order selection is to implement the identification algorithm 
recursively in the order. For scalar systems this method is often referred to as the Levinson-Durbin 
algorithm (LDA). By extending the conventional LDA to multidimensional systems, it can be used 
in the present invention resulting in a significant saving in computations compared to standard 
techniques. The recursion is stopped when the residuals are small enough, thus minimizing the 
number of necessary operations. 
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The present invention is computationally more complex than a conventional MLSE receiver, i.e. a 
receiver wherein an equalization method is used. If it is necessary to reduce the computational 
complexity, pre-filtering is used. For illustration, M=K+L is considered and a general ARX model 
is used. The idea is to concentrate the main energy in B into fewer taps i.e. reduce the time 
dispersion of the filter and then to use the model in the equalizer. The first row of equation 15 is 
then expressed as equation 21, where the introduced notation has obvious definitions and q'* is the 
unit delay operator. A pre filter F(q**) is introduced and the filter of equation 2*1 and equation 22 is 
obtained. The remaining rows of eq 15 are similarly transformed. The choice of the nimiber of 
remaining significant taps depends on a trade-off between how much we need to reduce the 
computational complexity and noise amplification. The determination methods are readily available 
in the signal processing literature. 

An embodiment of the invention is to combine the model of equation 15 with other receiver 
structures to reduce the computational complexity, while keeping the loss in performance at a 
modest level. An example is the DFSE algorithm earlier mentioned. The DFSE is well described in 
the literature. And is directly applicable to equation 15 by one skilled in the art. 

The invention can be used with one or more antennas. In equation 1 h would then be an mxL+1 
matrix where m is the nvunber of antennas. Correspondingly, in equation 4, h would be an 2mxL+l 
matrix and x(t) would have 2m elements. 

It is possible to get some additional performance if oversampling is introduced. Oversampling in 
this case refers to sampling faster than the bit rate. Typically the oversampling factor, denoted by P, 
would be equal to 2 or 4, but other values are of course possible. The extension to handle the case 
with oversampling is very similar to the extension to multiple antennas decribed above. When 
oversampling with a factor P, h in equation 1 would be a P x L+1 matrix. Row k of equation 1 
would correspond to samples k, k+P, k+2P,.... Correspondingly, in equation 4, h would be a 
2PxL+l matrix and the vector x(t) would have 2P elements. By introducing oversampling we get 
more channels than with symbol rate sampling, which gives larger potential for the interference 
rejection. On the other hand, the higher the oversampling factor, the higher the correlation between 
adjacent channels (in h). Therefore there is a limit to what performance gain can be achieved with 
oversampling. 
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In one embodiment of the invention prior knowledge of the receiver/transmitter and pulse shaping 
filter is taken into consideration. In general, when the modulation is linear, the channel is the 
convolution between the physical channel and a combination of the receiver, transmitter and pulse 
shaping filters. This latter combined filters is known to the receiver and it is possible to use this 
information in the demodulation process. It has been shown that the composite channel vector can 
be written as in equation 23, where G is a matrix formed from the known filters and c is a vector 
with the parameters of the physical channel. This formulation is particularly useful when we use the 
oversampling described above. The reason is that the different channels all consist of the same 
physical channels, but differ only in the matrix G which can be computed in advance. Hence we 
have greatly reduced the number of parameters in the problem. Receivers using multiple antennas 
can also benefit fi-om this kind of modeling. In this case the large dimensionality of the problem 
makes it important to keep the number of parameters to a minimum. 

Decision direction is a well-known method to improve the performance of a receiver in a 
communication system. The concept is based on the idea that once we have estimated all (or a 
certain nimiber) of the transmitted symbols, we make a hard decision about these symbols and use 
the decisions to re-estimate the model parameters with increased accuracy. This methodology, 
which can be repeated a desired number of times, is well documented in the literature. In 
combination with the method of the invention, decision direction is particularly iiseful. The reason 
is that in most practical communication systems, the number of training symbols is small in order 
not to waste information bandwidth. In the present invention, there are more parameters to estimate 
than in a conventional receiver, hence the quality of the parameter estimates can be significantly 
increased if decision direction is used. 

The preferred embodiement of the invention is described in Figure 2. The possibly over-sampled 
complex-valued signal , received with one or more antennas, is split into its real and imaginary part 
to obtain the vector in Equation 4. This vector is fed to the sequence estimator which determines 
the transmitted symbols using estimates of the parameters in the filters and of the model in 
Equations 15-18. These estimates are obtained using an available training sequence, and the number 
of parameters in the filters are obtained as an output of the order estimation algorithm. The model 
order can be obtained in several ways. One method described in the present invention relies on a 
test (described in 1-5 above) on the residuals from the parameter estimation, another method is 
decision direction from the metric obtained as an output of the sequence estimator; the model order 
giving the smallest metric when applied to the training sequence is used to estimate the data. 
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M,HD(0 + e(0 
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